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D , Abstract 

^ ■ 

We propose a method using perturbation theory in the running couphng 

constant and the idea of scahng to determine improved actions for lattice field 
K> ' theories combining Wilson's renormalization group with Symanzik's improve- 

j^ ' ment program . The method is based on the analysis of a single renormaliza- 

tion group transformation. We test it on the hierarchical 0{N) invariant a 
model in two dimensions. 



1 Introduction 

Improved actions for lattice field theories were proposed by Symanzik |S[ to reduce 
systematic errors due to a finite lattice spacing. The improvement parameters are 
calculated by bare perturbation theory. Improved actions also naturally arise in 
the block spin renormalization group approach of Wilson |[WK|| . There they come 
in form of a renormalized trajectory of lattice actions completely free of lattice 
artefacts. One way to approach the renormalized trajectory starting from any rea- 
sonable ansatz is to iterate renormalization group transformations. A different way 
due to Hasenfratz and Niedermayer |[HN|| for asymptotically free theories consists of 



a classical treatment of a single step. It gives the asymptotic form of the renormal- 
ized trajectory called perfect action. We propose to approximate systematically the 
renormalized trajectory in the spirit of Symanzik by a scaling analysis of a single 
renormalization group transformation in terms of perturbation theory in the running 
coupling constant. 

Consider the renormalization group flow of an asymptotically free lattice fleld 
on a unit lattice with action 5(0) generated by a renormalization group transforma- 
tion R. Imagine a one dimensional renormalized trajectory S^'^{(f), /) parametrized 
by the value of a running coupling /. A block spin transformation will do nothing 
to an action on this renormalized trajectory but change the value of the running 
coupling / into /'. The whole dynamics of the renormalization group is then en- 
coded in this change of the running coupling. (/ can for instance be marginally 
relevant, which means that f = f + A/?(/) with A(3{f) = (32p + 0{p) where 
(32 > 0.) An action on the renormalized trajectory therefore satisfles the discrete 
Callan-Symanzik type equation RiS'^^(0, /') = S^^{(l),f) and is said to scale. Re- 
call that we formulate everything in terms of transformations on unit lattices. (The 
reader who prefers to see explicit dependence on lattice spacings is invited to rewrite 
the equations using dimensional analysis.) 

In the domain where the running coupling is small, perturbation theory can be 
applied to compute S^'^{(j), f) (at least in a small fleld region) as a solution to the 
Callan-Symanzik type renormalization group equation. The strategy is as follows. 
Consider a lattice action of some general form 5(0) = J2aPa{f)Oa{(p). Here Pa(/) 
denote polynomials in the running coupling containing free improvement parameters. 
In the renormalization group terminology the action can contain irrelevant terms. 
We then compute a block spin transformation to a given order s of perturbation 
theory in the running coupling / to obtain RiS(0, /) = RiS*^'*-'(0, /) -|- 0{f'^~^^). From 
the result we determine the effective running coupling /' as a function of /, invert 
this relation, and express the effective action in terms of /'. We can then compare 
the original action as a function of / with the effective action as a function of /'. The 



idea is then to demand that the action be invariant up to corrections of order s + 1 to 
determine the form of the action and fix the values of the improvement parameters. 
Imagine that we have succeded in doing so. The result is the exact renormalized 
trajectory S^^\(l), f) up to order s in the running coupling. In other words we have 
removed all scaling violations to a given order in the running coupling. 

In practice the improvement procedure is performed step by step in perturba- 
tion theory. To get started we can take an unimproved action and see what kind 
of terms are generated in a first or second order calculation. Then the action is ex- 
tended by the corresponding terms with couplings depending appropriately on the 
running coupling with undetermined coefficients. Eventually this procedure has to 
be iterated until no further terms appear. Fortunately a finite order calculation can 
only produce a finite number of terms. Then the coefficients are fixed as described 
above completing the first improvement step. 

In a pilot study we have worked out this improvement program for the hierar- 
chical two dimensional 0(A^)-invariant cr-model. The existence of a renormalized 
trajectory in this model and asymptotic freedom (for A^ > 1) have been proved rig- 
orously by Gawedzki and Kupiainen ||GK|| . We have no doubts that our action can 
be proved to be an excellent approximation to the true renormalized trajectory in a 
small field region using [pK|| kind of bounds . As expected the leading term of the 



renormalized trajectory can be obtained by a perfect action calculation a la ||HN 
We take this perfect action as the starting point of the improvement program. The 
hierarchical cr-model itself is of marginal interest. However, it is an optimal labo- 
ratory to test renormalization group ideas. The implementation of our program to 
the full 0(A^)-invariant cr-model is under investigation. 

Let us mention that the improvement also works when dealing with more than 
one coupling to be renormalized. The renormalization group is then reduced to 
a transformation of a set of running couplings, the Callan-Symanzik type equation 
being exactly the one above. It is conceivable that the improvement program admits 
a nonperturbative formulation in terms of polymer representations. 

2 The hierarchical 0{N) model 

In the hierarchical (local approximation) renormalization group for general N com- 
ponent models we are left to deal with functions Z{(f)) of a single variable (p with 
values in IR . For 0{N) models one assumes further that Z(M(p) = Z{(j)) for all 
M G 0{N), requiring Z{(j)) to be a function of |0|. We consider the hierarchical 
renormalization group transformation in two dimensions given by (the case L = \pl 



of ^ ^, 

RZ(^)=Ar|d/i,(C)Z(^ + C)' 
It preserves 0(A^)-invariance. 7 is a real positive number and 

d/.,(C) = (27r7)^d^Cexp(^-^) (2) 

the corresponding Gaussian measure on H . A/" is an optional normalization con- 
stant. Concerning background material on hierarchical models see [ pK|| and refer- 
ences therein. In the following we will write 

Z(0)=exp(-r(0)). (3) 

A sufficiently general form of V^(0) for a models close to the renormalized trajectory 
proves to be 

v{<P) = Y.PaU)m-)T (4) 

a>2 J 

with P2(/) = 4"+0(/^) and PaW) = 0{f^) for a > 2. The running coupling is given 
by the inverse radius. In this setup is not restricted to take values on the sphere 
with radius j, the reason being that our recursion does not preserve this condition. 

The normalization constant M is conveniently chosen such that ^(0) = for |(/)| = j. 
We adopt this renormalization condition. When computing single renormalization 
group transformations we will speak of the previous action as the bare action and 
of the outcome as the effective or renormalized action. 



3 Perfect action 

To calculate the perfect action [|HN|] for the hierarchical 0{N) model we write the 



potential in the form 

V{(P) = r2V(^ - 1) (5) 



r 



with r = J. The hierarchical transformation then takes the form 
exp(-r'M^I^|-l))=Ar(^)"^/d^eexp('-^-2rMl^^ + e|-i: 



(6) 



The fluctuation field has been rescaled by r. It turns out in perturbation theory 
that the effective radius is 

r' = r-|(iV-l)i + 0(l). (7) 

In the limit when r and hence also r' is sent to infinity the fiuctuation integral can 
be evaluated by the saddlepoint method. Rescaling also the block spin ip hj r the 
equation for the perfect action becomes 

V^^m - 1) = ^inf^ l^-^e + 2V^^(|^ + el - 1)) . (8) 

This equation can be solved by the ansatz C2{\4>\ — r)^ with a single parameter C2- 
The solution is C2 = j-. The perfect a( 
trajectory in this model is then given by 



The solution is C2 = j-. The perfect action approximation for the renormalized 



V^^^(0) = f(-I0|-If, (9) 

47 r 

The right way to think of this formula is as a line of fixed points parametrized by 
r of the classical renormalization group transformation. Note that unlike ||HN|| the 



action is not just multiplied by r^. We have tested this approximation numerically 
as will be explained below. 

4 Perturbation theory 

The perturbation expansion for the 0{N) model can be computed to high orders 
using computer algebra. Let us explain the method in a second order calculation 
for the perfect action. As bare potential we take 

v{<j>) = ^{-\<p\-ir. (10) 

The effective potential will be 0(A^)-invariant. Without loss of generality we can 
therefore take the block spin to be given hj ip = {r + \l')e with e an iV component 
unit vector, say (0, . . . , 0, 1)^. The shift of \& serves to place us into the minimum 
of the bare potential. We then decompose orthogonally the fiuctuation field into 
C = (je + TT with respect to the direction of ip. The one component variable a is 
the radial fiuctuation field while the A^ — 1 component variable n is the tangential 



fluctuation field. The bare potential is expanded in powers of the coupling f = -. 
Up to second order it is given by 

47 27 47 47 47 

Z!^V - ^^an^ - ^aV + ^{nY + 0(f). (11) 
47 27 27 167^ ^ ^•^ ^ ^ ' 

One immediately observes a linear and a quadratic term in a which cannot be treated 
as perturbations. The mass term changes the a covariance to ^7. The source term is 
taken into account by a shift of the radial variable cr = ^ — |\t'. The renormalization 
group transformation then becomes 

exp(-V'(vI/)) = exp ('-i-vI/2W'|d/ii^(0 /d^,(7r) exp (-2V(v[/,e, vr)) (12) 

in terms of C, and vr. The bare potential takes the form 

V(^,e,vr) = i-M/TT^ + -^e^2 _ 

87 47 

^'.^V - ^^^vr' - — eV + -^{iry + Oif). (13) 



I67 47 47 I67 

At this point perturbation theory is applicable. Although the potential is non- 
polynomial to begin with only finitely many terms show up at finite order with a 
leading trilinear vertex. Note that the perfect action is recovered when fiuctuations 
are completely neglected. The perturbation expansion is straight forward using the 
Gaussian correlations 

/d/xi^(0e'"=(|)"n(2m+l) (14) 

m=0 

and 

. n—l 

J dfi,in)in'r = 7" n (2^ + ^ - !)• (15) 

m=0 

Computing the fiuctuation integral to second order perturbation theory we obtain 
an effective potential of the form 

V'(^) = (^ - ^(^ - 1)/') ^' + \(N - l)f^ + 0{f). (16) 



We then determine the value 5r of \& at which the effective potential attains its 
minimum and substitute \E' = $ + 5r. The change of r is due to the linear term in 
^ . The meaning of this variable is $ = j-?/'! — r' with r' = r + 5r the renormalized 
radius. To second order perturbation theory in / the change of the radius is 5r = 
— ^{N — 1)/ + 0{p). From this we find a renormalized coupling /' = ^ of the form 

f = f + l{N-l)f + 0{f)- (17) 

(The vanishing of the f^ term of the A/5 function follows from a fourth order calcu- 
lation.) In particular we have confirmed that the model is perturbatively asymptoti- 
cally free for A^ > 1. That is, when perturbation theory applies we find a flow where 
the coupling slowly grows and the radius slowly shrinks. The effective potential 
becomes 

V'($) = (^ - Yq{N - I) A $' + 0{f) (18) 

in terms of $ = j-?/'! — r' and /'. We also see that this action remains invariant in 
the sense of scaling to flrst order. Scaling violation shows up in a second order flow 
of the overall prefactor. In a zeroth improvement step they can be compensated for 
by changing the bare action into 

V'W=(^ + 4V)(H-i)'^ (19) 

The correct value of the improvement parameter is q = |(A^ ^1)- The resulting 
action can then be seen to scale even to second order. 

To appreciate this result consider for instance the polynomial bare potential 

n0) = ^(l0^-l)'. (20) 

Perturbation theory reveals that its effective action contains among other terms a 
cubic one at flrst order in the effective coupling. Therefore this action does not scale 
even to flrst order. 



5 Improved action 

Let us now also remove the scaling violations of second order. As a bare potential 
we take the second order improved one which we write in the form 

V^'\<P) = P?\f)(\<P\--\ (21) 



withPf (/) 



J- + |(A^ — 1)/^ . The effective potential computed to third order 
perturbation theory is given by 



V'(2)($) 



-XN-l)f 



12 



$2 + l(iV-l)/'3$3 + o(^M) ^22) 



in terms of /'. What is new is a cubic term in $ = |0| — r' . As it will be generated 
anyway to third order it is natural to include it to this order in the bare action. In 
other words let us make the ansatz 



y(3)(0) 



p?\^)m-j\ +i^f(/)(i0i 



/. 



(23) 



Mt 



>{2), 



.(3) 



,(3) 



(3), 



with Pr{f) = P2 if) + 4 7^ and Prif) = 4'{N - l)/'^. Here we have antici- 
pated a possible cubic correction to P2{f) which in fact turns out not to exist. The 
ansatz contains two new improvement parameters. Computing again the effective 
potential to third order, the correct values of these are seen to be C2 = and 
C3 = ^. This action reproduces itself up to fourth order corrections. The improve- 
ment scheme can now be iterated. Suppose that we have found the potential 



V^W(0) = ^PW(/) 



a=2 



f, 



(24) 



which scales up to order s. It's image under a renormalization group transformation 
evaluated in perturbation theory to order s + 1 is of the form 



\/(^+i)(0) = f:pi^+i)(/) 



a=2 



/. 



(25) 



containing polynomials Pi*'*'^-'(/) = Pa^Kf) + c^'^^^^*^^ for a < s and Pj+il/) 



c 



(s+l) f/s+1 



s+1 



f . This general form again is reproduced to this order. The corresponding 



effective potential to order s + 1 is given by 



s+l 



Y'(s+i)($) = ^ p<s+i) (/')$« + 0{f 



is+2\ 



(26) 



a=2 



With p^(^+^)(/) = [pW(/') + c;(^+i)r+i 



for a < s and P^^i 



(/) 



J(S + 1) fls + l 

-s+l J 



The effective coefficients c'j^^^^^ depend linearly on their bare counterparts c^'^^^\ 
(To order s + l they have no other choice.) In order that there be no scaling 
violation to order s + l the polynomials P'J^^^^\x) and Pjf^^\x) have to be equal. 



From this we obtain a system of linear equations for the coefficients c^*^^\ This 
system turns out to have a unique solution. To order /^ it is given by 



\2^ 



2 VJ^ 47 8^ V 56 12 ; •^ 

/ 57l7^(A^-l)^ 3271677^ (iV-1)^ 10305 7^ (A^ - 1) \ . 
V 84 ^ 70560 ^ 3472 J ■' 



_ 2'- 

36 '' '" \ 360 ' 336 ^ 

2^ 



pf(/) ^ ^(A--i)/'- P^^'J!-n^^^^':-^M f 



^ ^•'^ 112^ '-' \ 1736 35280 j '' ' 

Plt\f) = -^(^-1)/'- (27) 

We observe that the power series for Pa{f) contains only even (odd) powers of / when 
a is even (odd). Furthermore we observe that the signs of the coefficients alternate. 
The complete series is not expected to converge due to instanton singularities. It 
would however be very interesting to apply the machinery of resummation methods 
to a high order approximation to the renormalized trajectory. What remains to be 
said to order /^ is the recursion formula 

r^j,l^^f-ll^_Af,o(fy (28) 

for the running coupling. For the sake of completeness we also include the recursion 

|(iV-l)i + g(3(iV-l)^ + iV-l)l 
J- (l8(iV - 1)^ - ll(iV - if + 32(iV - 1)) - 
for the effective radius. 



6 Numerical results 

Fortunately we are not left alone with perturbation theory. Hierarchical renormal- 
ization group transformations can also be computed numerically. Therefore we are 



able to test our perturbative results and find the limits of their validity. Going back 
to the point of departure we write the transformation in the form 



exp(-V'(|^|))=Ar|d^eexpf-i-e^-2V(|^ + e|)V 



(30) 



with V(|0|) = V{(j)). The fluctuation integral is most conveniently done in N dimen- 
sional polar coordinates given by |^| = R and ip ■ C, = \ip\ R cos 6'. The independent 
polar angles are immediately integrated over, changing just the (irrelevant) normal- 
ization. The result is 

OO TV 

exp {-V{\^\)) =J\f' I dR j de w{R, e)B U\^l)\^ + 2\%l)\R cos 6 + rA (31) 



with a weight function 

w{R, e) := i?^-i sin^-2 9 exp ( -— ^' ) (32) 

and a Boltzmann factor (squared by the block volume two) 

5((p):=exp(-2V(^)). (33) 

^p denotes a positive real valued single component field (the modulus of ip). In a 
numerical evaluation the range of the i?-integration can be restricted to [0, Rmax] 
with Rmax sufficiently large. The contribution of {Rmax, cxd) is both supressed by 
the exponential decay of the weight function w{R,6) and the decay of the Boltz- 
mann factor. Furthermore it is sufficient to follow the flow of B{ip) in an interval 
f € [^min, fmax] with numerically negligible corrections. Let us emphasize that 
the interval should include a large field region not covered by perturbation theory 
in order to test our improved potential. We parametrize \/{f) by a set of N^ val- 
ues Vj = V((^j) at equidistant points ipi G [(pmin,fmax]- Intermediate values are 
computed by cubic-spline interpolation. The integral (0) is performed by stan- 
dard numerical methods yielding a set of effective values V^ from which we compute 
parameters r' and C'^ by a least square fit using the ansatz 

^'i^) = tc'aiv-rr. (34) 

a=2 

The renormalized trajectory can be computed by iterative renormalization group 
steps. (All numerical results refer to parameter values 7 = 1, A^ = 3, fmin = 

9 



r — 30, (fmax = r + 30, N^ = 160.) We start with an unimproved bare poten- 
tial V^'-°''(|0|) = 1/(47)(|0| — rY and iterate renormalization group transformations. 
After each step the resuhing effective potential comes closer to the renormalized 
trajectory giving successivly better approximations (see figure [|). We took the val- 
ues of the couplings C^(/') =: P^'^if) after 10 iterations as estimates of the true 
renormalized couplings at the radius r' = 1//'. By varying the bare radius r we 
computed the renormalized trajectory for a set of renormalized radii r' G [0, 100]. 
The results were then compared with the coefficients P^\f') of the improved per- 
fect action approximation V'^^\ As expected one finds excellent agreement (fig. ^) 
for large r = ^. The relative deviation of the perturbative coefficients from their 
nonperturbative values is less than 10~^ for r > 30 and remains smaller than two 
percent for r > 5. The improvement compared to the unimproved perfect action is 
especially visible in the large field region (fig. |). 

Around r = 5 the deviation increases dramatically. This breakdown of pertur- 
bation theory in the crossover region from the ultraviolet to the infrared regime is 
illustrated in fig. p. The renormalized trajectory connects the UV-fixpoint with the 
IR-fixpoint. The IR-fixpoint potential is a quadratic single well whereas the UV- 
fixpoint potential has a double well shape. Along the renormalized trajectory the 
radius r and the depth of this well shrink simultanously. At some (pseudo) critical 
point Tc ~ 2.04 (see fig. ^ the radius and the depth vanish and remain zero all the 
way to the IR fixed point. Obviously the coupling / = - diverges at this point and 
looses its meaning as expansion parameter. Though highly interesting in its own 
we do not intend to investigate the infrared behavior in terms of the fiow in this 
paper. Let us only mention that it is conceivable that one could performe a change 
of coordinates to another expansion parameter in the crossover region. A natural 
candidate is the mass parameter in the potential since the double well becomes very 
flat before it turns over. 

Acknowledgment 

We would like to thank M. GrieBl for helpful comments. 

References 

[HN] P. Hasenfratz and F. Niedermayer, Perfect lattice action for asymptotically 
free theories, preprint BUTP-93/17 

[GK] K. Gawedzki and A. Kupiainen, Continuum limit of the hierarchical 0{N) 
non-linear a-model, Comm. Math. Phys. 106, 533-550 (1986) 

10 



[S] K. Symanzik, Continuum limit and improved action in lattice theories, (I) 

Principles and 0^ theory, Nucl. Phys. B226 (1983) 187-204; (II) 0{N) non- 
linear a model in perturbation theory, Nucl. Phys. B226 (1983) 205-227 

[WK] K. Wilson and J. Kogut, The renormalization group and the e expansion, 
Phys. Rep. C12, No. 2 (1974) 75-200 



0.25010 
0.25000 
0.24990 
0.24980 
0.24970 



0.24960 

0.02000 




0.02008 



Figure 1: Flow of the coupling C2{f) towards it's renormalized value 
P^'^if) under consecutive applications of the RG transformation. The 
calculation was performed with a bare radius of 50. The other cou- 
plings Ca{f) with 2 < a < 6 behave similar. 
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Figure 2: Comparison of the perfect renormalized coupling -P^^(/) 
with it's improved and unimproved approximations P2 (/),-P2 (/) 
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Figure 3: Large field behavior of the effective potential. 
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Figure 4: Flow of the effective potential. The plot shows V(|-?/'|) after 
n=0,3,6,9 and 30 RG steps starting with the unimproved perfect action 
approximation \/*^°^(|-?/'|) = l/(47)(|-?/'| — r)^ at r = 5 and ending in the 

HT fixpoint \/^t(|^|) = 1/(47)|^P. 
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Figure 5: Flow of the radius in the JR. At Vc ~ 2.04 the effective 
radius vanishes with r' ^ Ir — rpl^'^. 
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